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Viable models of modified gravity designed to produce cosmic acceleration at the current epoch, closely mimic 
the ACDM model at the level of background cosmology. However, this degeneracy is generically broken at the 
level of linear perturbations, where the modifications induce a peculiar scale-dependent pattern. A similar pattern 
is expected in models of coupled dark energy. 1 present the main results published in Pogosian and Silvestri [I], 
on the growth of structure in f(R) theories of gravity, and in Zhao et al. 2\, on the potential of upcoming and 
future tomographic surveys to detect departures from the growth of cosmic structure expected within General 
Relativity with a cosmological constant. 



1. Introduction 

Observations strongly favor a universe that has 
recently entered a phase of accelerated expan- 
sion [314] , This poses a puzzle for Modern Cos- 
mology since standard General Relativity (GR), 
applied to a universe which contains only radi- 
ation and dust, has difficulties fitting the data. 
The cosmological constant A, a dark energy (DE) 
fluid (i.e. a yet unknown component with a nega- 
tive equation of state) as well as modifications of 
GR on large scales are the different proposals to 
address the phenomenon of cosmic acceleration. 

It has become increasingly evident that tests 
based solely on the expansion history cannot 
place significant bounds on many models that ad- 
dress cosmic acceleration. Indeed, although cur- 
rent geometrical probes place tight constraints on 
the expansion history, models of modified gravity 
in general have enough freedom to reproduce any 
desired expansion history. In other words, at the 
level of background cosmology there is a degen- 
eracy among different models of cosmic accelera- 
tion [1171819110111] . However, this degeneracy is 
typically broken at the level of cosmological struc- 
ture formation; indeed, models of modified grav- 
ity that closely mimic the cosmological constant 
at the background level can still give significantly 
different predictions for the growth of structure. 
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An instructive example is offered by f{R) theo- 
ries of gravity. These models are degenerate with 
ACDM (i.e. GR with cold dark matter and a 
cosmological constant) at the level of background 
cosmology. However, they predict a significantly 
different dynamics of perturbations, with a pe- 
culiar scale-dependent pattern for the growth of 
structure. Studying in detail the formation of 
structure in this class of models allows one to gain 
insight on the growth of perturbations in gener- 
alized models in which a scalar field, (being it the 
dark energy field or the additional degree of free- 
dom introduced by modifications of gravity), is 
coupled to dark matter. 

From this analysis it becomes clear that the 
large scale structure of the universe offers a 
promising testing ground for GR. It is there- 
fore important to explore to what extent one 
can detect departures from GR in the growth of 
structure with present and upcoming cosmologi- 
cal data. In this regard, I review the results of 
Zhao et al. [2], where we have used parametric 
forms to quantify the potential time- and scale- 
dependent modifications of the growth of struc- 
ture. Since the effects of the modifications of 
gravity on the growth of structure are typically 
equivalent to, or can be mimicked by, a more ex- 
otic dark energy field, (i.e. one that couples to 
dark matter and/or carries anisotropic stress), we 
have focused on the distinction between ACDM 
and any alternative. 
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2. f(R) theories of gravity 



scalaron 



In 1979 Starobinsky showed that a de Sitter 
phase in the early universe can be achieved by 
adding to the Einstein action a function f(R) oc 
R 2 |12ll3j . More recently, f(R) theories have 
been revisited in the context of cosmic acceler- 
ation, starting with 14 15J and followed by many 
others. The underlying idea is to add functions of 
the Ricci scalar which become important at late 
times, for small values of the curvature. Let us 
consider the action for a generic f(R) model 

S = i / <*V=5 [£ + /(«)] 

+ d 4 x^g£ m [x t ,g^}, (1) 



where x% represents the matter fields, which 
are minimally coupled and, therefore, fall along 
geodesies of the metric g^y. The field equa- 
tions obtained varying the action (JTJ) with respect 
to g^ v , and the energy- momentum conservation 
equations are, respectively, 
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where /r = df/dR. The Einstein equations @ 
are now fourth order in the scale factor. The addi- 
tional dynamics introduced by the modifications 
of gravity can be seen explicitly by taking the 
trace of Einstein equations, to get 



□/« = I (R + 2/ - Rfu) 



k2 ( 1^- dv ° s 



(4) 



which is a second order equation for with a 
canonical kinetic term and an effective potential 
V c s(fii). Therefore the modifications of the Ein- 
stein action introduce an additional scalar degree 
of freedom and the role of the scalar field is played 
by fit, dubbed the scalaron [13]. One can asso- 
ciate the following Compton wavelength to the 
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The scalaron adds extra dynamics to the theory 
and mediates a fifth-force between matter fields. 
Both these features can be problematic for the vi- 
ability of the f(R) model. Consequently, special 
care needs to be taken in order to ensure a stable 
high-curvature regime, to have a positive effec- 
tive Newton's constant and to satisfy local tests 
of gravity |8ll6ll7ll8ll9l20l2lj . After all the nec- 
essary conditions are imposed, viable f(R) mod- 
els closely mimic ACDM at the level of expansion 
history, with undistinguishable departures from 
w e s = — 1. Nevertheless the predictions for the 
growth of structure can differ significantly from 
those of ACDM, as we will see in the following 
subsection. 

2.1. Dynamics of Linear Perturbations 

Let us concentrate on linear scalar perturba- 
tions to a flat Friedmann-Robertson- Walker uni- 
verse, in Newtonian gauge. In particular, since 
we are interested in the growth of structure, let 
us focus on cold dark matter and consider linear 
sub-horizon scales. In Fourier space, the equa- 
tions governing the evolution of perturbations are 
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where a prime indicates derivation w.r.t. In a, 6 = 
Sp/p is the density contrast of dark matter and ^ 
and $ are, respectively, the perturbations to the 
time-time and space diagonal metric component. 

Comparing eq. ([7]) with , it is clear that 
the modifications have introduced a Yukawa-type 
correction to the Newtonian potential, with char- 
acteristic lengthscale A c ([5]): the scalaron me- 
diates a fifth force. Therefore, the Compton 
wavelength of the scalaron introduces a scale 



3 



k ( hi Mpc) 



0.5- 
0.2- 
0.1- 
0.05- 
0.02- 
0.01- 




4 


















1 1 1 I 11 1 I T 1 



I I I I I II I I I I I 

54321 05 43210 



z z 



0.8 0.9 



Figure 1. The evolution of the growth factor 
for the CDM [A(k,a)/a]/[A(k,ai)/ai\ as a func- 
tion of redshift and scale. The left panel corre- 
sponds to an f(R) model with w B s — — 1 an d 
A° re 10 2 Mpc, the right one to ACDM. The 
dashed line crossing diagonally on the left plot 
corresponds to the Compton scale (0 ■ 



which separates two regimes of sub-horizon grav- 
itational dynamics during which gravity behaves 
differently, as can be noticed in Fig. [TJ On scales 
A A c , the scalaron is massive and the fifth force 
is exponentially suppressed, thus deviations from 
GR are negligible. However, on scales inside the 
Compton radius, the scalaron is light and devi- 
ations are significant. The relations between <E> 
and 'J, and the relation between them and the 
matter density contrast, will be different below 
the Compton scale and that affects the growth 
rate of structures. In particular, on scales be- 
low A c , the modifications introduce a slip between 
these potentials, leading to \& ~ 2<F For modes 
that cross below A c during matter domination the 
effect of modifications is maximized as the sum 
of the potentials grows in the absence of back- 
ground acceleration. When, however, the back- 
ground starts accelerating, the potentials begin 
to decay but at a lesser rate than in the ACDM 
model. A characteristic signature of f(R) theo- 
ries would be a non-zero Integrated Sachs- Wolfe 
effect (ISW) during the matter era. 

Overall, we observe a scale-dependent growth 
pattern. The modifications introduced by f(R) 
models are analogous to those introduced by more 
general scalar-tensor theories and models of dark 
energy coupled to dark matter. The dynamics 



of perturbations is richer, with different observ- 
ables described by different functions, not by a 
single growth factor. Therefore combining differ- 
ent measurements, such as galaxy counts, weak 
lensing and their cross-correlation with the ISW 
effect, we can test gravity on cosmological scales. 
But what is the power, of upcoming and future 
tomographic surveys to constrain modifications 
from ACDM? In [2J, the authors addressed this 
question employing a parametrization to describe 
the modifications, i.e. the rescaling of the New- 
ton's constant and the gravitational slip. In what 
follows I review the main results of [2j. 

3. Searching for modified growth patterns 

The linearized Einstein equations provide two 
independent equations relating the metric poten- 
tials and matter perturbations: the Poisson and 
anisotropy equations. Just like in f(R), eqs. (J7J 
and ([8]). in models of modified gravity and more 
exotic models of dark energy, the relation between 
the two Newtonian potentials, and between the 
potentials and matter perturbations, can differ 
from the ACDM case [22I23I24| . We parametrize 
the changes to the Poisson and the anisotropy 
equations as follows: 

k 2 * = -^rv(a,k) P A (9) 

§=7(a,*0, (10) 

where pA = pS + 3 3 ^-(p + P)v is the comoving 
density perturbation, p,(a, k) and 7(0, k) are two 
time- and scale-dependent functions encoding the 
modifications of gravity and/or the contribution 
of an exotic dark energy fluid. In ACDM p, = 1 = 
7- 

To arrive at a suitable parametrization of the 
functions p(a, k) and 7(0, fc), we note that mod- 
els of modified gravity typically introduce a tran- 
sition scale which separates regimes where grav- 
ity behaves differently. For instance in f(R) the- 
ories the functions p and 7 are equal to unity 
at early times and on large scales (i.e. on scales 
that are larger than the Compton scale, A c , of the 
model), and they transition to a modified value 
on smaller scales and late times. This time- and 
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scale-dependent transition can be mimicked via 
the following functions 



(i(a,k) = 



l+/3iAf k 2 a s 

1 + \\ k 2 a s 
1 + /3 2 A| k 2 a s 



1 



X 2 k 2 a s 



(11) 
(12) 



where the parameters Af have dimensions of 
length squared, while the pi represent dimen- 
sionless couplings and s > 0. Expressions (fTTj) 
and (p~2|) coincide with the scale-dependent 
paramctrization introduced in [25 . It is easy to 
see that this parametrization corresponds to the 
addition of a Yukawa-type interaction between 
dark matter particles, with an effective coupling 
— 1) and a range Ai. In the case of f(R) indeed 
we have seen that the modifications introduce a 
fifth-force mediated by the scalaron. 

In order to determine how well upcoming and 
future surveys can constrain the Modified Growth 
(MG) parameters {s, fit, 02, \x, X2} we employ 
the standard Fisher matrix technique [29] . 

As our observables, we use all possible two- 
point correlation functions (both auto- and 
cross-correlations) between the Galaxy Counts 
(GC), Weak Lensing shear (WL), and Cos- 
mic Microwave Background (CMB) temperature 
anisotropy, across multiple redshift bins, in ad- 
dition to the CMB E-mode polarization autocor- 
relation, and the CMB E-mode and temperature 
cross-correlation. 

In order to compute the theoretical expecta- 
tions for these observables, we assume fiducial 
values for the parameters and evolve the pertur- 
bations using MGCAME0 g, (a modified version 
of the publicly available CAMB [26127] ). For the 
fiducial values of the parameters we consider sev- 
eral models based on f(R) and Chameleon theo- 
ries. All f(R) models correspond to a fixed cou- 
pling, Pi = 4/3 (= 8/3 p 2 ) and to s ~ 4. The only 
free parameter is the lengthscale A2 — y/JhXi. 
Here I report results for the specific choice of 
A| = 10 4 Mpc 2 . In addition to the effects of mod- 
ified gravity, the evolution of all cosmological per- 
turbations depends on the background expansion. 
We restrict ourselves to background histories con- 
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sistent with the flat ACDM model, since that is 
currently the best fit to available data and since 
popular models of modified gravity give an ex- 
pansion history which is effectively the same as 
in ACDM; the main differences between models 
indeed arise from the evolution of cosmic struc- 
ture. We do not, however, fix the values of the 
cosmological parameters Q m h 2 , flbh 2 and h, nor 
the spectral index or the optical depth. We vary 
them along with the modified growth parameters, 
using the WMAP 5-year data best fit [28] for their 
fiducial values. 

Finally, we assume CMB T and E data from the 
Planck satellite [32], the galaxy catalogues and 
WL data by the Dark Energy Survey (DES) [3U] 
and Large Synoptic Survey Telescope (LSST) 31j, 
complemented by a futuristic SNe data set 
provided by the ongoing Canada-France-Hawaii 
Telescope's Supernovae Legacy Survey [33|34j . 
the Nearby Supernovae Factory [35136] , and a fu- 
ture Joined Dark Energy Mission space mission, 
such as the Supernovae/ Acceleration Probe [37] . 
We describe the specifics of the experiments in [2] . 

3.1. Results 

We evaluate the Fisher errors on the MG pa- 
rameters, as well as the usual set of standard cos- 
mological parameters. Fig. 13.11 shows the 68% 
C.L. contours for the fiducial case of an f(R) 
model with Compton scale A 2 , = 10 4 Mpc 2 . This 
is just one representative case among the sev- 
eral fiducial models considered in our analysis [5] ■ 
In general, we find that the five parameters are 
correlated and the degeneracy between Ai and 
A2 is strongly positive, while Pi and Pi have 
a strong negative correlation. This is expected 
since one can enhance the growth by either raising 
Pi, which increases the effective Newton's con- 
stant, or lowering /J 2 , which enhances the rela- 
tive strength of the Newtonian potential which 
drives the clustering of matter. This also explains 
the degeneracy between the two length scale pa- 
rameters Ai and A2. 

Comparing the constraints on the MG param- 
eters from different datasets, we see that the 
constraining powers of WL, GC and WLxGC 
are comparable and that they provide much 
more stringent constraints than WLxCMB and 
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Figure 2. The 68% confidence contours for the 
five MG parameters for the fiducial f(R) model 
with A| = 10 4 Mpc 2 , as constrained by several dif- 
ferent combinations of correlation functions from 
various experiments (shown in the legend). Each 
of the ellipses is plotted after marginalizing over 
all other parameters that were varied. 



GCxCMB. This is, in part, because there are 
many more different correlations between WL 
and GC than there are cross-correlations with 
CMB. Also, correlations of GC and WL with 
CMB suffer from a larger statistical uncertainty, 
since only the ISW part of the total CMB 
anisotropy is correlated with large scale structure. 

Comparing the constraints from DES with 
LSST, we find that, typically, LSST can improve 
constraints on individual MG parameters by a 
factor of 3 or better. With all tomographic data 
combined (no priors used), DES give relative er- 
rors (i.e. the la marginalized errors divided by 
the corresponding fiducial values) on the MG pa- 



rameters of order 20%-40%, while LSS can go be- 
low 10% levels. 

4. Conclusions 

Models of modified gravity and more exotic 
dark energy typically induce a scale-dependent 
growth pattern and a time- and scale-dependent 
slip between the gravitational potentials. Study- 
ing in detail the class of f(R) theories, we have 
learned how the modifications come into play and 
which peculiar signatures they imprint on the LSS 
and CMB. In general, the dynamics of pertur- 
bations is richer, breaking the degeneracy that 
characterizes cosmic acceleration models at the 
background level. Weak lensing, galaxy count, 
the ISW effect and their cross-correlations offer a 
powerful testing ground for GR on large scales. 
It is therefore important to determine to which 
extent current and future data can detect depar- 
tures from GR in the dynamics of perturbations. 

In [2] we have used a five-parameters descrip- 
tion for the rescaling of the Newton constant, 
/z(a, k), and for the ratio of the metric poten- 
tials, 7(0, k), to investigate the power of upcom- 
ing and future tomographic surveys to constrain 
departures from GR. From the point of view of 
scalar-tensor theories, e.g. f(R), these parame- 
ters are related to the coupling in the dark sector 
and to the characteristic mass scale of the model. 
We have then studied in detail the constraints on 
the five-parameters based on several fiducial mod- 
els, expected from Weak Lensing, Galaxy Counts, 
CMB and their cross-correlation spectra as seen 
by Planck, DES and LSST (additionally using 
Planck and SNAP to constrain the standard set 
of cosmological parameters). We have found that 
DES can provide 20%-40% level constraints on in- 
dividual parameters, and that LSST can improve 
on that, constraining the MG parameters to bet- 
ter than 10% level. Overall, with DES, and espe- 
cially with LSST, one will be able to significantly 
reduce the volume of the allowed parameter space 
in scalar-tensor type models. 

From this analysis we have learned that up- 
coming and future surveys can place non-trivial 
bounds on modifications of the evolution of per- 
turbations, even in the most conservative case, i.e. 
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considering only linear scales as we did. These re- 
sults are model-dependent, but they motivated us 
to pursue model-independent methods such as the 
Principal Component Analysis [38 ; we recently 
reported the first results in [39 . 
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